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OUTLINE 


• Motivation 

• PDF modeling of reactive flows 

• The Lewis PDF module 

• Validations and applications 

• Current research 

• Technology transfer 


COMPUTATION OF TURBULENT COMBUSTION 
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GOVERNING EQUATIONS 


P.t + {P u i),i = 0 

(pUi).t + (P u j u >).] = ~P.i + r, ;j 

{pE) %t + ( pUjE)j = + ‘I* 

(pYk),t + (pujYk)j = (pDYkj)j + u k 


- 4 , 



Aj 


- dA 

dxj 


CLOSURE PROBLEM: 


u"u" 


Yyu'j 


P‘Wi 


U{ — Ui + tt;, 

y. — Y + Y' 

i l — i i . J- i ) 


— Turbulence Modeling 

Analogy of shear stress: Jifiusion model. 

??? 


pwi = pw(Y\, ...,Y n ,T ) 

But in generai: .. 

jjWi # fni>(Y\, y„,T ) 
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PDF Modeling of Turbulent Reactive Flows 
Current status 

« Assumed PDF (Spalding, 1971; Gosman &: Lockwook. 
1973; ...) 

0 Advantage: simple, fast. 

0 Disadvantages: Need unique mixture fraction; assumed 
shape may not be real. 

« Composition PDF (Pope, 1976; Dopazo &c O’Brian, 
1974) 

o Advantage: Reaction rate treated exactly; existing mo- 
ment closure codes easily adapted. 

0 Disadvantages: Turbulent diffusion needs model. 

t __ 

. Velocity-Composition joint PDF (Pope V Chen 19S0, 
Pope 1981) 

o Advantage: Reaction rate treated exactly; no d Fusion 
model needed. 

0 Disadvantages: Models for velocity field relatively un- 
tried; Require more computer resource. 


PDF Modeling of Turbulent Reactive Flows 

• Objective: 

0 Develop models that can accurately simulate finite rate 
chemical reactions in turbulent nows. 

0 Develop and validate independent PDF modules. 

0 Technology transfer. 

* Criteria 

0 Accuracy and robustness. 

0 Practical in terms of today’s computing power. 

0 Easy integration with existing industry computational 
platform. 
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PDF Modeling of Turbulent Reactive Flows 


• Approach: 

0 Joint pdf method for scalar compositions. 

0 Moment closure schemes for velocity field. 

0 Develop hybrid solver consisting of Monte Carlo method 
and finite-difFerence/linite- volume method. 


PDF Modeling of Turbulent Reactive Flows 


. Current status (Lewis) 

(pP),t + (P < Uj\Yi, li > P)j -I -{pwjP)' Yj 


o Continuous mixing model developed. 

0 Model for compressibility effect proposed. 

0 2D and 3D Monte Carlo PDF module developed. 
0 Validation studies. 

0 Code released l.o industry during a workshop. 
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Validation Cases 

• Scalar field in homogeneous turbulence. 
. Oblique shock. 

• 2D supersonic hydrogen combustor. 

• Axisymmetric supersonic combustor. 

« Piloted flame near extinction. 
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Scalar field in homogenous turbulence 
pdf compared with Gaussian distribution 


Current model 


Modified curl's model 




u:c 2. AiyiOptoi*« pdf distribution for #. iCAi*r in lionw 
turbu!««. — preacai rjoiicl; - • * C*u«.xn. 


t\j yiftpioiit pdf a.om (of * icaIa* in hotno- 

* V :;i30‘,‘ i lurbuicMCc ■— mocifi'cd Cu-! r"<5<ie!. - * * C/4U*iiAii 


Scalar field in homogenous turbulence 
3rd and 4th moments compared with 
Gaussian 


Current model 


Modified curl's model 
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Evolutiou of moaienl* L-cm the present ootid. — 
it^odjird deviation, • * ■ O.i x fourth ccaifii moocal, - - 
0.0 1 x tilth tcitiril rr.oineoi, #0.1 x fourth moment for 
C*uS9tu distribution. Q 3.01 x lix'.a moment for Cius#nn 
distribution . 


involution of inO:i«c:itj from the modified Cud 
tTtodd. — Jtiadifd dcv:ii:oa, • ■ - 0.01 -x :oiir*.b erntui 
••trace;:’.. - - 0.0001 x lixtii ecotrii triomeat. o 0.01 x fourth 
i:iui:ic 3 t far Causiiau distribution. C 0 0001 x sixth niorneiil 
for Cauuijui distribution. 
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Temperature across an oblique shock: 
pdf solution compared with analytical 
prediction. 
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Supersonic hydrogen combustor 
(Exp. Burrows & Kurkov, 1973) 



Temperature Contour (pdf solution) 
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Supersonic hydrogen combustor 
Mole fraction: 

pdf solution compared with exp. data 



Coaxial burner: geometry and test condition 
(F.xp. Cheng, et al. 1991) 
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Mean H20 mole fraction 
Coaxial burner 

pdf solution compared with exp. data 



Piloted flame (Masri et al. : 1994) 

Fuel: 45% CO, 15% H2, and 40% N2 
Flame close to extinction 



255 


Piloted Flame 
Mean Temperature 



Piloted Flame 
Flame Location 
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Current Projects 


0 Application of PDF module to emission predictions 
0 Incorporate general chemistry procedure. 

0 Incorporate spray models. 

0 Use parallel computing for the PDF module. 


Collaboration with industry and 
technology transfer 

. Features of independent pdf module: 

0 Easily coupled with any existing industry flow codes, 
o Novel averaging scheme to reduce memory requiement. 
o General chemistry package. 

0 Parallelized workstation version. 

. Technology transfer: workshops 

0 July, 1993; code released to 15 US institutions. 

0 October. 1994. 
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